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Abstract 

Following the linear programming prescription of Ref. [6] , the d <8> d Bell diagonal 
entanglement witnesses are provided. By using Jamiolkowski isomorphism, it is shown 
that the corresponding positive maps are the generalized qudit Choi maps. Also by ma- 
nipulating particular d(&d Bell diagonal separable states and constructing corresponding 
bound entangled states, it is shown that thus obtained d <8> d BDEW's (consequently qu- 
dit Choi maps) are non-decomposable in certain range of their parameters. Keywords: 
Entanglement witness, Bell state, Generalized Choi map. Non decomposable 

PACs Index: 03.65.Ud 
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Entanglement is one of the most fascinating features of quantum mechanics. As Einstein, 
Podolsky and Rosen [1] pointed out, the quantum states of two physically separated systems 
that interacted in the past can defy our intuitions about the outcome of local measurements. 
Moreover, it has recently been recognized that entanglement is a very important resource in 
quantum information processing [2]. A bipartite mixed state is said to be separable [3] (not 
entangled) if considered as a convex combination of pure product states. 

One of the approaches to distinguish separable states from entangled ones involves the so 
called entanglement witness (EW) [4] . An EW for a given entangled state p is an observable 
W expectation value of which is nonnegative on any separable state, but strictly negative on 
an entangled state p. 

There is a correspondence relating entanglement witnesses to linear positive (but not com- 
pletely positive) maps from the operators on Hilbert space H A to the operators on Hilbert 
space Hb via Jamiolkowski isomorphism, or vice versa [5]. 

Here in this paper, by using the prescription of Ref[6], i.e., reducing the manipulation of 
generic Bell-state diagonal-entanglement witnesses to an optimization problem, we find d® d 
Bell states diagonal entanglement witnesses (BDEW). By using Jamiolkowski isomorphism, we 
show that corresponding positive map is the generalized qudit Choi map. Also by manipulating 
particular qudit Bell diagonal separable states and constructing corresponding bound entangled 
states, we show that thus obtained d® d, BDEW's (consequently qudit Choi maps) are non- 
decomposable in certain range of their parameters. 

The paper is organized as follows: 
In section 2, we introduce the the d®d Bell diagonal separable states and corresponding PPT 
states. In section 3 we give a brief review of entanglement witness and we show that finding 
generic Bell states diagonal entanglement witnesses for d®d systems reduces to a LP problem. 
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In section 4 and 5 we find generalized d <g> d Choi entanglement witnesses of first and second 
types, respectively. Section 6, is devoted to investigation of non-decomposability of generalized 
Choi entanglement witnesss. The paper is ended with a brief conclusion. 

2 Bell diagonal d<S) d separable states 

Here in this section we introduce some sets of the d(E)d Bell diagonal separable states(BDSS). 
In general the set of BDSS consists of the following three categories: 
Set 1: The first set of BDSS can be constructed from Bell states 

yd l=0 

where u = and k = 0, l,...,d— 1, simply by summing over the following Bell states 
projection operators: 

El^oX^ol = \ E d5 hV \l)\l) ® (2-2) 

fe=o a 1,1' =0 

where we obtain one of the BDSS of first type, defined as 

Po = E |VW(^o| = E 10(^1 ® 10(^1- (2-3) 

fc=0 1=0 

Acting Shift operator on (2-3) we can find the other BDSS of first type as follows 

p m = (I d ® S rn ) E |^o)(^o|(/d ® S m )\ m = 1,2,.., d-1. (2-4) 

k=0 

Therefore, the first set of BDSS consists of the following d BDSS 

d-l d-l 

p m = ^\^km){^km\ = ^\l){l\®\l + rn){l + m\, m = 0, 1, .., d - 1. (2-5) 

k=0 1=0 

Set 2: 

In order to obtain the second set of BDSS, we need to consider the following sum of 
projection operators 

€ = E k) E fa) ® E(^l E (v m \, (2-6) 

i=0 j=0 1=0 m=0 
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Vd 

Now summing over the free index k and doing some routine calculations we obtain the 
following BDSS of second type 



d-l d-l d-l j d-l ^ 



Po = E & = E(E -rjW + k)) ® (E -ri(j\(j + k\), (2-7) 

fc=0 k=0 i=0 V « j=0 V« 

where it can be written in terms of Bell states as 

Po^El^XV'ofel- (2-8) 

fc=0 

The remaining BDSS of second type can be obtained by applying powers of Modulation oper- 
ator over (2-7) separable state, 

P'm = E( fim ® WofcX^K^™ ® ^ m = 1, 2, d-l. (2-9) 

fc=0 

Hence, the second set of BDSS consists of the following d BDSS 

d-l d-l 

p' m = El^)(Vw| = E w m( - i - i) \i)\i + k)®{j\(j + k\, m = 0,l,..,d-l. (2-10) 

fc=0 i,j,k=0 

Set 3: Finally in order to obtain the third set of BDSS, we first define the following set of 

vectors 

\u i ) = a i (u k ) i \i), (2-11) 
\w i ) = b i (u k ) i \i), (2-12) 
then using above vectors, we construct the following sum of projection operators 

€ = (E^(^)i^>)(E^(^) J '0'l) ® (X>(^I0)(E MO m M), (2-i3) 

i=0 j=0 i=0 m=0 

where we assume that a, = 6j and a i+1 = us~ m (ii. Now, summing over free index k and using 
the relation) | 2 = and doing some straightforward calculations we obtain the following 
BDSS's of third type 

pl= £ co nki ^=\i)\i + k)®^=(j\(j + k\)u;-^, (2-14) 
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where it can be written in terms of Bell states as 

d-i 



Pn = EM^IWfelM^' (2-15) 
k=0 

Therefore, the third set of separable states are defined as 

d-l 
k=0 

where p% = p' . So far we have introduced (3d-l) separable states. Other separable states 
can be obtained by applying powers of local Shift (S) or Modulation(f2) operators on these 
separable states. For instance we consider 3 <8> 3 system which consists of 8 separable states. 
Now, we can obtain some other separable states simply by acting powers of Shift operator on 
above defined separable states as 

(J 3 <g> S 2 ))p'i{h ® S 2 Y = |10)(10| + |02)(02| + |21)(21|, 
(h ® S*))p>i(h ® = |20)(20| + |02)(02| + |11)(11|, 
(J 3 ® S))p'l(I 3 ® Sy = |01)(01| + |12)(12| + |20)(20|, 
(J 3 ® SM(/ 3 ® 5) f = |01)(01| + |22)(22| + |10)(10|. 
It is easy to see that some of thus obtained BDSS's lie at the boundary of separable states 
[6]. One can show that the p l s = po,p| = Po = Po are orthogonal to the optimal W re d, i.e., we 
have Tr[W re dPs} 1,2. Hence, p l s ,i = 1,2 lie at the boundary of the separable region. 

Also convex sum of these states i.e, p\ = Xpg + (1 — A)p| is orthogonal to the optimal W re d, 
i.e., we have Tr[W red p\] = 0. Hence, p\ lies at the boundary of the separable region. 

At the end of this section we try to introduce some positive partial transpose (PPT) opera- 
tors states which will be used later, in the investigation of non-decomposability of generalized 
Choi entanglement witnesss. We consider the following d <g> d density matrix defined as 

1 — p 

PPPT = P(|V ; 00)(V ; 00|) H rf - ^ 1 ? 1 + ^ 2 P 2 + ••• + P'd-lpd-l), (2-18) 

where it is positive for < p, /xi, p,2, ■, ■, Pd-i < 1 and Yltzl P-% — 1- Now in order to make the 
partial transpose of density matrix (2-18) to be positive, i.e., to obtain PPT density matrix, we 
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need to use the relations pJ A = pi and substitute the partial transpose of Bell state projection 
operator (|^oo) (^ool) 

a m,l=Q 

in partial transpose of Eq.(2-18). Now, the positivity of partial transpose of density matrix 
(2-18) implies that 



p T / > Q ^ p < max i min i_tj_ || fc / i / - = 1 (2-20) 

For yUi = p 2 = ... = yUd-i = ^rr the parameter p is optimal and equal to p — ^, where d® d 
density matrix reduces to 

l_ pd-l y 
PPPT = WIV'OoXV'Ool) + — ^2f*iPi],P< ~ d - ( 2_21 ) 

As we will show later in section 6, for certain range of parameter p the density matrix (2-18) 

become entangled (actually bound entangled due to its PPT property). 

In remaining part of this section we consider particular case of d = 3 

1 — p 

Pppt = pQVW^ool) H g— (A»iPi + P2P2), (2-22) 

where the positivity of its partial transpose implies that 



p<max{ ^Elt^ }i (2-23) 

where for p\ = p 2 = \ the parameter p is optimal and it is equal to p = |. 

Finally we construct 3 <8> 3 PPT state of second type by using (2-10) separable states as 

1 — p 

Pppt = p( I ^00)^00 1) + — ^ — (MiPi + P2P2), (2-24) 

where it is positive for < p, p±, p 2 < 1 and Yh=\ Pi — 1- 

Now, the positivity of partial transpose of density matrix (2-24) implies that 

T 2p\ - 2 + 2p\ - 2p x p 2 + <$P2 + 3/ii + 3 J 4 + Qpip 2 - < &p\ - 3/x| 

p p PT > =>- p < max\ — 7, 7, — ~r — : }. 

hppt- v- x 2 {8 + p 2 2 - p ± p 2 + pj + 3p ± + 3p 2 ) 1 

(2-25) 

For p\ = p 2 = \ the parameter p is maximum and equal to p = \. 
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3 Bell state diagonal d® d entanglement witness 

In this section we give a brief outline of some of the main features of Bell state diagonal d <S> d 
entanglement witness together with linear programming prescription for finding it. For further 
information reader is referred to Ref. [6]. 

Let S be a convex compact set in a finite dimensional Banach space. Let p be a point in 
the space with p which is not in S. Then there exists a hyperplane that separates p from S 
[5, 7, 8]. 

A hermitian operator (an observable) W is called an entanglement witness (EW) iff 

3p such that Tr(pW) < (3-26) 

Vp' e S Tr(p'W) > 0. (3-27) 

Using these definitions, we can restate the consequences of the Hahn-Banach theorem [8] 
in several ways: 

1- p is entangled iff there exists a witness W such that Tr(pW) < 0. 

2- p is a PPT entangled state iff there exists a non-decomposable entanglement witness W 
such that Tr(pW) < 0. 

3- a is separable iff for all EW Tr(Wa) > 0. 

From theoretical point of view this theorem is quite powerful. However, it is not useful to 
construct witnesses that detect a given state p. 

We know that a strong relation was developed between entanglement witnesses and positive 
maps[5, 7]. Notice that an entanglement witness only gives one condition (namely Tr(W p) < 0) 
while for the map (I a <8> 4>)p to be positive definite, there are many conditions that have to be 
satisfied. Thus the map is much stronger, while the witnesses are much weaker in detecting 
entanglement. It is shown that this concept is able to provide a more detailed classification of 
entangled states. 
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Following Ref[6], one can expand any trace class observable in the Bell basis as 

W= £ W iii2 \i> iii2 )(4> iii2 \ (3-28) 



where \ip ili2 ) for (0 < i\ < d, < i 2 < d) stands for the orthonormal states for a d®d Bell state. 
After some calculations similar to those of [6], the trace-1 Bell state diagonal W observable 
can be written as 

W = r^ + (l-r) J2 q iii2 \^ iii2 )^ iii2 \. (3-29) 

»1»2=° 

The observable given by (3-29) is not a positive operator and can not be an EW provided 
that its expectation value on any pure product state is positive. For a given product state 
1 7) = I a) 1 1 a) 2 the non negativity of 



implies that 



7V(W|7)(7|) > (3-30) 



d 2 y^ 1 a P 

Tl IT < r < 0, (3-31) 



1 - d 2 E , q P 

^i 1 i 2 =0 ^»1*2 H 

where P. —\< 7 I 1b. . >\ 2 . 

As for the completeness of the Bell state E ili2 \i J i l i 2 ){i'i 1 i 2 I = 1> the determination of r c 
reduces to the following optimization problem[10] 

minimize C~ = T\. . q P (7) 

0<P ili2 (7)<J ( 3 - 32 ) 
£. ■ P 4 4 (7) = l- 

n j 2 l i l 2 v ' ' 

Always the distribution P . satisfies < P (7) < -, for all pure product states[6l. One can 
calculate the distributions t (7), consistent with the aforementioned optimization problem, 
from the information about the boundary of feasible region. To achieve the feasible region we 
obtain the extreme points corresponding to the product distributions P 4 s (7) for every given 
product states by applying the special conditions on q ili2 S parameters. C 7 themselves are 
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functions of the product distributions, and they are in turn are functions of 7. They are not 
real variables of 7 but the product states will be multiplicative. If this feasible region constructs 
a polygon by itself, the corresponding boundary points of the convex hull will minimize exactly 
C 7 in Eq. (3-32). This problem is called LP , and the simplex method is the easiest way of 
solving it. If the feasible region is not a polygon, with the help of tangent planes in this region 
at points which are determined either analytically or numerically one can define new convex 
hull which is a polygon and has encircled the feasible region. The points on the boundary 
of the polygon can approximately determine the minimum value C 7 from Eq.(3-32). Thus 
approximated value is obtained via LP. 



4 Generalized qudit Choi map of first type 

Compared with qutrit Choi positive map [9] (f)(a,b,c) : M 3 — > M 3 , Generalized qudit Choi 
map of first type <f)(a , ■ ■ ■ , a d _i) : M d — > M d is defined as 



'ao,---,ad- 



,(P> 



^ a p n + a 1 p 22 H h a d „ x p dd ... 

Od-iPn + a p 22 H h a d - 2 pdd ■■■ 



y ... aipn + a 2 p 2 2 H h a p dd 

-p (4-33) 

where p G M d . Using Jamiolkowski [5] isomorphism between the positive map and the opera- 
tors, we obtain the following d <g> d entanglement witnesses corresponding to Choi map 

Wchoi = 

Y d-l d-l d-l 

— ■ ■ — (oo lV , fco)(V ; fco|+ai l^fciX^fciH ^ a d-i l^fc,rf-i)(V'fc,rf-i|-^l^oo)(^oo|), 

d{a H h a d _i - L) k=Q k=Q fc=0 

(4-34) 



T / \ T 

• 2 .. u d - 1 
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where it can written in terms of separable states (2-5) in the following form 

Y d-l 

Wchoi = -j, ■ ■ 7t(E a mPm - d\ i/jqo) (ipoo I ) • (4-35) 

rf(a + • • • + a d _x - 1) ^ 

For product state (7) = |a)i|a) 2 = 73 ( 1 u uj 2 ... u/ _1 ) ^73 ( 1 u u 
the non negativity of (4-35) implies that 

«o + oi + ... + a d -i > d. (4-36) 

Similar to BDEW we expand | - 0oo)('0oo| using the identity operator and the other Bell diagonal 
states: 

d-l 

|-0oo> (^00 1 = I# - E l^vKV'yl- (4-37) 
Then we reduce EW to the following form 

j d-l 

W^ftoi = -77 ; ; 7r(-(d- a )I d 2 + d V |^fco)(^fco| 

a(a + h a rf — ij fc=1 

d-l d— 1 

+(ai + d- oq) l^fciKVfcil H h (ad-i + d- ao) E IV'm-iKV'm-iD- ( 4 "38) 

fc=0 k=0 

Comparing with BDEW (3-29) we have 

r = d ( d - °°) (4-39) 

(ao + --- + a d _i-l)' 

and the EW operator is defined as 

1 2 (1 — r) d_1 d l 

Wchoi = 7T— 7; ~T\ — 7 7 7 l^o)(^fco| + (ai+rf-a ) ^ |^ fc i><Vfci| 

d 2 {{d 2 - 1) + (1 - d)a + ai H ha«j_i) ^ ^ 

d— 1 

+ • • • + (a d _! + d - ao) E iV'M-iXV'M-il)- ( 4 "40) 

fc=0 

Note that if r is negative, as introduced in EW above, this operator will be positive, but not 
a completely positive map. Using (4-36) inequality, the minimum negative eigenvalue of choi 
EW (4-40) is given by 

^ + (l-r)— — — ■ ■ ■ -<0 , o mjn = mm{a ,ai,-,fli-i}) 

rf 2 d((rf 2 - 1) + (1 - d)a + ai H ha d _i) 

(4-41) 
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where above inequality is satisfied for r < and 1 < a < d. 

By using (3-30) for non-negativity of the observable W c h i we find the distributions as a 
function of qij. The minimum value of C 7 is obtained from the boundary of the feasible region, 
i.e., we have 

1 



(C 7 ) = 



{{d 2 - 1) + (1 - d)a + ai + • • • + a d _i) 



(p 1+ (ai+ '- a °W - • •+ (arf - i+ /~ ao) ^), 



ri 



ri 



(4-42) 



where Pi = X^=i P k0 ^2 = Y,t=o P k i anc ^ *Pd — Y,i=l P k d-i- We can nn d the extreme value of 
(Pi , • • • , Pd) which is obtained under the product states (7) = |a)i|a) 2 as 

Pi =| a | 2 | fa | 2 + I ai | 2 | A | 2 + ••• + I a d _i | 2 | A*-i | 2 



-3 II a || A) I + I «i II A I e^ 1 + •••+! || A-i | | 2 
P2 =| « | 2 | A | 2 + I ai | 2 | A | 2 + • • • + I ad-i | 2 | A) I 2 



Pd =| «o | 2 | A-i I 2 + I "i | 2 | A) I 2 + ••• + I "d-i | 2 | A-2 | 2 



(4-43) 



where = 



/ \ 



y otd-i ) 



and I a) 2 = 



y A-i y 



. One can obtain the extreme points of the 
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{V u V 2 ,---,V d ) as 



( 1 ^ 



l«)l = l«)2=73 



« 1 



a; 

2 



CO 



V w ) 



ah = 



\ / 



v / 



V / 



where u = e <* . 



, |«>2 = 



V / 
' o ^ 





lit h 


v ; 



v / 



(^l = i^2 = i--,^=J) 



(Pi = 0,P 2 = 1,P 3 = p 4 = • • • = v d = 0) 



(Vi=V 2 = --- = Vi-i =0,Vi = l,V i+ i = ■ ■ ■ = V d = 0) 



(4-44) 



The convex combination of all extreme points provide a convex or a feasible region, then 
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we have the following optimization problem 



minimize (C 7 ) 
subject to 



(((2 2 -l)+(l-ci)a„+ai+-+a d _ 1 ) 1 ci 

1 d-1 ' 1 d-1 ' 2 ' 3 



(ai+d-ao) p^ _|_..._|_ ( a d-i+ rf ~ a o) p^^ 



1 - 3=1^1 - ^T^3 - ^2 + • • • + V d < 

1-V 1 -V 2 V d > 

V u V2,---,V d >0. 



(4-45) 



Analytically, we have been able to show that we will have violation only from the hyperplanes 

d 

(d-1)- dVi -(d- l)Vi ~(d-2) ^ Vj = , i = 2, ...d. 

Now let us assume that the maximum value of the violation from the planes is A < 1. Thus, 
the equation of the plane passing through the new extreme points, parallel to the above plane, 
is obtained. Next we derive the intersection of the following adjacent planes 

dV x + (d- 2)V 2 + (d - 1)7> 3 + • • • + (d - l)V d - (d - 1 + A) = 



dVx + (d- 1)V 2 + (d-l)V 3 --- + (d- 2)V d - (d - 1 + A) = 
V x + V 2 + • • • + V d - 1 = 



Vi = 
Vi = l 



, i — 1, d 



(4-46) 



, i — 2, d 
Vi+V 2 + --- + V d = ^ 

where new extreme points are obtained from intersecting the above hyperplanes. Next we 
calculate C 7 for all the newly obtained extreme points and compare them with one another. 
Some easy calculations give the minimum value of the parameter C 7 which is independent 
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from A: 

d — \ 1 -|~ a m»n-«0 

d {{d 2 - 1) + (1 - d)a + ai H h a d _i) 

(4-47) 

then the critical value of the parameter r is obtained as 

d C-v 



l-d 2 C^ . 

frrnn 



(4-48) 



For a = ai = • • • = a^-i = 1 the parameter r reduces to r c = —d corresponding to the well 
known reduction map 

W red = J2 Pi -d\i>J(i> 00 \. (4-49) 

i=0 

On the other hand, EW (4-40) must have positive trace under any product state |7)(7|. Thus 
the introduced r in (4-39) must satisfy 

did — a ) — d 2 (X, . , , 

r>r c ^- 1 0) > Z,™" , 4-50 

(ao + • • • + a d _i - 1) 1 - cPC 7mm 

where the inequality is satisfied for all value of 1 < a < d — 1. 

5 Generalized d<g)(2 Choi entanglement witness of second 
type 

Similar to (4-35) one can define the second generalized d <g> d Choi entanglement witness by 
second set of Bell diagonal product states p' m as 

d-l 

W choi = a iPi - MJi^oo I • ( 5 " 51 ) 

Comparing with (4-38) entanglement witness W' Choi reduced to 

Y d-l d-l 

W Chd = 7T— I ~ a °^ 2 +d T, \^ok)(^ok\ + + d - a )p' k ). (5-52) 

d{a H h a d _i - 1) k=1 k=1 
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Comparing with BDEW (3-29) we have 

d(d-ao) (5 _ 53) 



(a H h a d _i - 1)' 

By using (3-30) for non-negativity of the observable W' choi we find the distributions as a 
function of g^. The minimum value of C'^ is obtained from the boundary of the feasible region, 
i.e., we have 

(r \ \ t^i , fai + d - op) ^ , , ( a d-i + <j ~ go) 

1 7j " ((d 2 -l) + (l-rf)a + a 1 + --. + a d _ 1 ) l/l+ " d 2+ ' ' ' + " d ^' 

(5-54) 

where 7^ = Efc=l ^ ^2 = Efc=o P i fc and = Efc=o p d -i, fc - We can find the extreme value of 
^P'li'P'n ' ' ' I'P'a) which is obtained under the product states (7) = |a)|/3) as 

V[ = |(E&I ^A +fc (Eto^^A + ^) - \ I Eto I on || Pi I I 2 , 00 = 

^2 = ^(Efc7=o a 'A+fc(Ef=o^ l ««+iA+fe+i) . . 

(5-55) 

n = ^(E£io^A +fc (Eto 1 ^ (d - 1)i «wA + ^) 
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where la) = 



/ \ 

a 



and \I3) 



\ O-d-1 J 



( 8 X 



. One can obtain the extreme points of the 



l«> = 1/3) 








V / 



a) 



i 

y/d 



,1/3) 



1 



V 1 / 

t 1 ^ 



V 1 / 



CO 



A-l 



•10 = 71 



i 



( i ^ 



V 1 / 



J 



1 



V 1 / 



d ' ' 2 



^ = 0) 



(K = o,^ = i,n 



7Vi = 0,7" d = l). 



(5-56) 



The convex combination of all extreme points provide a convex or a feasible region which is 
completely similar to feasible region of Wchoi entanglement witness (4-34) and then optimiza- 
tion problem is similar to previous optimization problem and finally C lmin is equal to (3-32) 
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and critical value for r parameter is equal to (4-48). 
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6 Non-decomposibility condition for generalized Choi 
entanglement witnesses 

As it is explained in first section, an entanglement witness W is non-decomposable iff there 
exists a bound entangled state p where Tr(pW) < 0. First we assume that Choi's EW is 
decomposable, hence it can be written as a convex sum of a positive operator P and partial 
transpose of a positive operator Q as follows 

W Chm = Q Ta + P- (6-57) 

To do so consider (4-49) reduction EW (W re d) which is a optimal decomposable EW (For more 
details see [6, 11]) hence it can be written as partial transpose of positive operator Q i.e., 
W red = Q Ta where Q = Ef=o pi - E^Lo ^l^)!^™,^)- Now usin S the reduction map one 
can decompose the Choi's EW of first type as 

W Choi = \W red + (l-\)P, (6-58) 

where positive operator P is 

~ . ciq — d A . | \ / / i / ^ \ v — i \ / i 

(E i= o a i - l ) d d (Ei=o a i - l ) d d ( d - X > k=l 

The above positive operators is Bell states diagonal, hence its eigenvalues, i.e., the coefficients 
of Bell states projection operators must be positive, therefore we should have 



d — a 
(Eto Oi - 1) 

aj( d ~ 1) 
(Eto 1 Oi - 1) 



A > 7^1-^ ( 6 - 60 ) 



A< ^Ti. \, , j = 0,...,d-l. (6-61) 



Generalized qudit Choi maps 



19 



Now, combining the inequalities (6-60) and (6-61) yields 

at > , i = l,...,d-l. (6-62) 

Also, by raining both sides of inequality (6-60) to (d-1) power and multiplying both sides of 
inequality (6-61), we get following inequality 

(d - a )^ < ^ < ^gjzig , ^ 



(Eto <k - ly- 1 ~ " " (Eto^-l)^ 1 
The above inequality implies that the Choi's EW is decomposable as long as its parameters 

satisfy the following inequality 



(d-g ) d - 
(d-1) 



aia 2 ...a d _i > — — j^zr- (6-64) 



For particular case of 3 <E> 3 Choi EW, Wchoi = apo + &Pi + C P2 — 3|'0 OO )('0 OO | above inequality 
reduces to 

bc > ^~ A a) \ (6-65) 

which is in agreement with [9]. 

If one of the inequalities (6-60) and (6-61) violates, then operator P will not be positive 
operator any more, consequently one cannot say that whether Choi EW is decomposable or 
non-decomposable, because decomposition (6-58) is not unique. In these cases, in order to 
find non-decomposability conditions for d <E> d Choi's EW (4-35), we try to detect d <g) d PPT 
density matrix (2-21) by them, namely we should have 

d-1 

Tr[W C hoipppT\ = (1 - p)Q2 OifM) + p(ao - d) < 0, (6-66) 



which yields the following lower bounds for parameter p 



. Ej=i 0>i[l>i /„ 

V > , ~ ^5=i • (6-67) 

Now, combining this lower bound with the corresponding upper one, due to PPT property of 
d® d density matrix (2-21), namely p < |, we get the following rang for parameter p 

Ei=l a ipi . 1 I r ro \ 

-T-. <P<- 7 , (6-68) 

d + E"=i a iPi - a o d 
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therefore, non-decomposability condition is 

d-l 

(d- 1)^2 ciiHi < d-dQ. (6-69) 

i=i 

First we consider 3 <E> 3 systems Choi EW where non-decomposability condition (6-69) reduces 
to 

bfn + cfi 2 <^^: (6-70) 



Now choosing ^ — 1 , fj, 2 — we get 



fe<V' (6 - n) 



on the other hand, the EW (4-36) condition implies that c > Similarly by choosing 

fj,i — 0,/i 2 — 1 together with (4-36) we get 

3 — a 3 — a . . 

b >— » c< ^-- ( 6 " 72 ) 

Summarizing above results, we can deduces that Choi's witness is decomposable if its param- 
eters satisfy (6-64) and it is non-decomposable otherwise. Since violation of the inequalities 
(6-60) and (6-61) will be equivalent to non-decomposability conditions provided that it remains 
an EW. Since these conditions can be the same as the non-decomposability (6-69) simply by 
appropriate choice of /x^, i = 1,2. 

Again one can conclude that d <E> d entanglement witness is decomposable if its parameters 
satisfy (6-64), other wise it is non-decomposable. Discussion about the decomposability or 
non-decomposability conditions of d <8> d Choi's EW is similar to 3 <8> 3 case. Obviously if 
the parameters Oj,z = l,...,d — 1 satisfy condition (6-64) it is decomposable but if some of 
the inequalities (6-60) and (6-61) violate, then one can show that, they are equivalent to 
non-decomposability condition. 

Decomposability condition of Choi's EW second type is similar to first type since, one write 

W Chm - xw red = ( (E ^~! 1)d + - d )\A )(A \ + ( (Stto i%_ 1)d " dp^i)) E l^oJ 
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d-i 



cii X 



(EtoV-l)d d(d-l) 




(6-73) 



i=i 



where the positivity of above operator yields results similar to (6-64). 

Discussion about its non-decomposability is hard in general. Since it is hard to find a PPT 
state consisting of separable state. So we restrict ourselves in the remaining part of this section 
only to 3 ® 3 EW second type. Now evaluating Tr\Wchoip'ppT\ with p' PPT given in (2-24), we 
get 



we see that non-decomposability condition is similar to (6-69), therefore decomposability and 
non-decomposability conditions of EW of second type is the same as with the first one. 

7 Conclusion 

It is shown that finding generalized d®d Choi Bell states diagonal entanglement witnesses can 
be reduced to an LP problem. A large group of non decomposable Choi entanglement witnesses 
have been defined by using Bell diagonal product states and corresponding d <g> d PPT states. 
We hope that in this way, one can study the Optimality and non-decomposibility of Choi 
entanglement witnesss for generic bipartite d\ ® di systems and multipartite d\ <S> c?2 <8> •■• <S> d n , 
which are under investigation. 
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